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Abstract 

By using multi-dimensional particle-in-cell simulation, we present a new regime of stable proton 
beam acceleration which takes place when a two-specie shaped foil is illuminated by a circularly 
polarized laser pulse. It is observed that the lighter protons are nearly-instantaneously separated 
from the heavier carbon ions due to the charge-to-mass ratio difference. The heavy-ions layer 
extensively expands in space and acts to buffer the proton layer from the Rayleigh- Taylor-like (RT) 
instability that would have otherwise degraded the proton beam acceleration. A simple three- 
interface model is formulated to qualitatively explain the stabilization of the light-ions acceleration. 
Due to the absence of the RT-like instability, the produced high quality mono-energetic proton 
bunch can be well maintained even after the laser-foil interaction concludes. 

PACS numbers: 52.40Nk, 52.35. Mw, 52.57.Jm, 52.65.Rr 
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In the past decades, plasma-based ion accelerators have attracted a 
to their potential applications for particle acceleration, medical therapy 



ot of attention due 
proton imaging, 
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and inertial confinement fusion (ICF) [2J. One of the most important issues is the devel- 
opment of laser-driven protons for radiation therapy of deep-seated tumours [3|. Numerous 
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experimental and theoretical studies have been devoted to producing such proton beams 
(|. However, their qualities such as beam collimation, energy spread (~ 20%), and peak 
energy (~ 58 MeV), are still unsatisfactory U 

Recently, with the rapid development of the laser technology, ultraintense ultrashort 
ultraclean (3U) laser pulse and ultrathin solid target have been extensively exploited to 
investigate the ion acceleration. One of the most straightforward acceleration mechanisms, 
radiation pressure acceleration (RPA) [6], is being revisited. The first RPA experiment [3] 
has shown that both the beam quality and the energy conversion efficiency are greatly 
improved. However, the ultrathin foil is very susceptible to the transverse instabilities [8], 
similar to Rayleigh- Taylor-like (RT) instability in ICF. It sets in at the very beginning of 
laser-foil interaction and develops at the unstable interface at the rate of a few laser cycles 
Gradually, the surface of the foil becomes corrugated by the laser radiation and the entire 
target is torn into many clumps and bubbles js)]. The final energy spectrum of the ions shows 
a quasi-exponential decay with sharp cut-off energy. Unlike the electron acceleration in the 
bubble regime jl0|, a stable proton beam acceleration in the realistic three-dimensional (3D) 
geometry has not been demonstrated either theoretically or experimentally. 

In this Letter, we report on a new regime of stable proton acceleration, where a two- 
specie ultra-thin shaped foil is illuminated by a circularly polarized (CP) laser pulse. We 
assume the heavier (lighter) ions to be carbons (protons), respectively. Particle-in-cell (PIC) 
simulations indicate that the RT instability only causes the spreading of the carbon ions. 
The protons, which are rapidly separated from the carbon ions, are "buffered" by the carbon 
ion cloud by riding on the stable proton-carbon interface. We demonstrate that, even though 
the RT-unstable carbon-vacuum interface is strongly deformed, the feed-through of the RT 
instability into the RT-stable proton-carbon interface is small. Due to the absence of the 
RT instability, the compact proton layer remains well-collimated even after the laser-foil 
interaction concludes. In order to elucidate the detailed acceleration process, we first describe 
the results of ID simulations. Discussion of the influence of the RT instability on the ion 
acceleration follows, backed up by 3D simulations. 



When a relativistic laser pulse illuminates a two-specie foil with thickness of a few wave- 
lengths, a collisionless shock wave is often excited and can efficiently accelerate the ions 
to high energies 



llfl. With the decreasing foil thickness, the laser radiation pressure com- 
petes with the shock wave and becomes strong enough to push the entire foil forward. As 
a result, the foil acceleration is dominated by the RPA. The critical foil thickness can be 



approximately estimated by 
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~ -—A, (1) 
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where a = eEi/m e cu is the dimensionless laser amplitude, m e the electron mass, n e the 
electron density, and n c is the critical plasma density c, u, A, and El are the light speed 
in vacuum, the laser frequency, wavelength and electric field, respectively. In the ID RPA 
model, the target motion equation is governed by 

P dt 27tc1 + /3' 1 ' 

where p = Y^- m^L is the target area mass density, and n, are the ion mass and density, 



(5 = v/c is the target velocity and 7 = 1/y 1 — /3 2 is the relativistic factor. We can see that 
the target dynamics is defined by the area density, not the detailed foil composition. In 
principle, the heavier ions can be efficiently accelerated to the same velocity as the lighter 
protons and electrons. 

I I 

We simulate the described mechanism using the PIC code VLPL [13|. The longitudinal 
length of the ID simulation box is x = 60A sampled by 6 x 10 4 cells, enough to resolve 
the expected density spikes. Each cell contains 100 numerical macro particles in the plasma 
region. The target is 0.1A thick, located at x — 10A and composed of carbon ions and protons 
with the same number density 46.7n c , which gives the electron density n e = 320n c . A CP 
laser pulse with the wavelength A = 1.06/im is incident on the target from the left boundary. 
The wave front of the laser arrives at the target surface at t = 10T , where T = A/c is the 
laser cycle. The laser pulse is homogeneous in space but has a trapezoidal profile (linear 
growth - plateau - linear decrease) in time. The duration is tl = 10To (ITq — 8Tq — ITq). 
The dimensionless laser amplitude a = 100 is chosen to satisfy Eq. ([TJ). 

Figure Ufa) shows the particle density distribution at t = 20T . In the initial stage, the 



laser pressure is transferred to the electrons, resulting in the charge separation 



12j. Because 



carbons and protons are initially co-located, the protons experience a higher acceleration 
due to their higher charge to mass ratio (Zi/rrii). The time for protons to separate from 
the carbon ions is approximately t sep = a/ 2Lm# / eEi = 2.5/s, which is so short that can 
be considered instantaneous. Later on, the two ion species start experiencing very different 
acceleration field, as shown by the red curve in Fig. [1(a). The considerably higher electric 
field inside the carbon layer compensates for carbon's lower Zi/rrii ratio enabling them to 
catch up with the protons. Eventually, both species travel together, without separating any 
further. The entire foil acceleration proceeds until the end of the laser-foil interaction at 
t = 35To. Fig. 1(b) exhibits the phase space distribution. We can see that the carbon ions 
fall back behind the protons, accompanied by a long low-density tail. The fact that both 
ions show an obvious "spiral structure" in phase space provides a direct evidence for 
the acceleration process described above. 
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Figure 1. (Color online), (a) Density distribution of the electrons (black), protons (green), and 
carbon ions (blue) at t = 20Tq. The red curve shows the electric field E x , which is normalized to 
Eq = m e cuj/e = 3.2 x 10 12 V/m. (b) The corresponding phase space distribution at t = 20Tq. (c) 
Energy spectrum of the carbon ions (blue, dark) and protons (green, light) at t = 30Tq. (d) Energy 
evolution in time from PIC simulations and the ID theory. The laser pulse is incident from the left 
side and touches the foil at t = WTq. 



The ion energy spectrum is shown in Fig. [TJ(c). At t — 30Tq, the peak energy of the 
carbon ions is up to 480MeV/u. For protons, all of them are accelerated to high energies 
although the energy spectrum is somewhat wider. Fig. [U(d) plots the ion energy evolution. 
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Here, we make use of the averaged energy for both species. At t = 35Tq, the laser-foil 
interaction is over so that the ion energy doesn't increase any more. Overall, the observed 
ion acceleration in the ID simulations is consistent with the predictions of the ID theory of 



Ref. 
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However, multi-dimensional simulations exhibit a radically different acceleration dynam- 
ics because multi-dimensional effects - such as transverse expansion of the bunch and the RT 
instability - come into play. In order to extend the ID model to the 2D simulations smoothly, 



we employ a shaped foil target (SFT) 15| to compensate for the transverse profile of the 
laser pulse. Taking the Gaussian laser for example, the foil thickness should be matched 
transversely by the Gaussian function L = max[L max exp(—y 2 joW L cut ], where L max is the 
maximal foil thickness, L cut the cutoff thickness, and Gt the spot radius. In the following 
2D case, the simulation box is X x Y = 80A x 32A, sampled by 16000 x 400 cells. Each cell 
contains 100 macro particles in the plasma region. The foil is initially located at x = 10A 
with parameters L max = 0.1A, L cut = 0.05A, and ax = 7A. The carbon ion density is 51.9n c , 
intermingled with protons of the density 8.64n c so that the total electron density is 320n c . 
A Gaussian laser pulse with the focal size <jl = 8A is incident from the left boundary. All 
other parameters are the same as in the ID case. 

Fig. |2]^a) shows the space distribution of the carbon ions and protons at different times. 
In each frame, the cyan color marks the carbon ions and the red color shows the protons. 
Obviously, the carbon ions behave totally different as compared with the ID simulations. 
They spread widely in space and do not form a compact mono-energetic bunch. On the 
contrary, the protons from the center part of the foil always ride on the carbon ion front and 
form a compact bunch. The sharp front separating the species is well defined and remains 
stable even after the laser-foil interaction concludes. We can get further understanding of 
the acceleration process from the phase space distribution in Fig. [3] (a) and (b). On the 
one hand, the carbon ions evolve into a wide cloud in space. On the other hand, their front 
co-moves with the fast protons so that the gap between the two species is always small. 
The protons show a clear spiral structure, like a "matchstick", which coincides with the ID 
simulation result in Fig. [U(b). We don't observe any obvious transverse instability in the 
compact proton layer. Fig. EJ^b) shows the proton energy spectrum. As expected, the peak 
is well pronounced and the dispersion is suppressed. The peak energy evolution is presented 
in Fig. []Jd), which is also in accordance with the ID RPA model. Fig. [3]^c, d) plots the ion 




Figure 2. (Color online), (a) Contours of protons (red, dark) and carbon ions (cyan, light) in the 
2D case at different time points: t = 30To,50To, and 70Tq. The colorbar represents the proton 
numbers log(N). (b) Proton energy spectrum at: t = 30Tq (square, black), t = 50Tq (circle, green), 
t = 70To (triangle, red), and t = 90To (star, blue). For comparison, Frame (c) shows the proton 
density distribution in a pure Hydrogen foil and Frame (d) corresponds to its energy spectrum 
evolution. 

energy-divergency distribution at t — 30T . The high quality proton bunch with the energy 
~ 500MeV and the opening angle ~ 5.5° forms and persists in time even after the laser-foil 
interaction concludes. 
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Figure 3. (Color online). Phase space of (a) protons and (b) carbon ions at t = 30To. An obvious 
spiral structure is observed only for protons. Frame (c) and (d) show the carbon ion and proton 
energy distribution as a function of the divergency angle at t = 30Tq . 
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The stability of the proton acceleration in the 2D simulations can be attributed to two 
effects. Firstly, the protons completely separate from the carbon ions and form a thin shell, 
which is a prerequisite for the stable proton acceleration. Such a separation of the ion species 
can be understood within the ID formalism developed in Ref. |l2| and the ID simulation 
above. Secondly, the heating of the carbon ions forms an extended cloud that prevents 
short-wavelength perturbations from feeding through into the thin proton shell. We can use 
a simple three-interface model as shown in Fig. H] to explain the stabilization. It is helpful 
to consider the problem from the classical RT instability 16| which occurs when a light 



fluid is accelerated into a heavy fluid. In the accelerating reference frame of the foil, the 
perturbation pressure p satisfies: 

d 2 

—Sp = -kl T 5p, (3) 

where /crt is the wavenumber of the RT-unstable mode. Noting that Sp is discontinuous 
across the unperturbed boundary, we obtain a solution 5p = Aie~ k ' RTZ + Bie k ' RTZ away from 
interfaces, with Aj and B t being the amplitude coefficients of the perturbation inside the 
layer consisting of the i'th species. In our case, both species have two interfaces: one with 
vacuum and one with the other specie. For carbon ions (i = C), the only unstable interface 
is the carbon-vacuum boundary, where the laser pulse interacts directly with the carbon 
plasma. We derive from this model that the amplitude of the RT instability is exponentially 
decaying away from the unstable interface: 

4^ ~ e - k ^ Lc , (4) 

where Lq is the thickness of the carbon ion layer. In the simulations, Lc is several times 
longer than Lh so that the long-wavelength perturbation in the carbon layer would take 
much more time to grow (recall that the the growth rate of the RT instability 7 oc ^/g/Xur, 
where g is the target's acceleration and A^r is the perturbation wavelength). The feed- 
through from the unstable carbon-vacuum interface to the proton layer is exponentially 
attenuated according to Eq. 03]). This simple qualitative argument explains the stability of 
the sharp carbon-proton interface. For the thin proton layer, it is also stable because the 
protons are much lighter than the carbon ions. Eventually, the entire proton layer is free 
from the RT instability. Besides, we believe that the small transverse size of the foil also 



benefits the stabilization of the proton acceleration in this case. Approximately, the minimal 
perturbation wavenumber can be estimated by /crt = 27r/y, where Y is the transverse foil 
size. At t — 30T , the carbon shell thickness is Lq ~ 10A and Y = 32A. Eq. already 
indicates a considerable suppression of the perturbation feed-through. 




Figure 4. (Color online). Schematic of laser-foil interaction in (a) C-H case and (b) pure H case. 
The red color marks the protons and the cyan represents the carbons. In case (a), there are three 
interfaces: carbon-vacuum, carbon-proton, and proton-vacuum. Only the first interface is unstable. 
In case (b), both interfaces are finally unstable. 

Now we compare the stable multi-component foil case with the pure proton foil case, 
where the RT instability is obvious. We again employ a matched SFT. All the parameters 
are the same as above except now n# = 320n c and the carbon ions are absent. Fig.^c) shows 
the proton density distribution in space. We can see that the foil disrupts gradually and two 
proton bunches with a lower density valley in the middle form. This is very characteristic 
for the RT instability driven by the laser radiation. According to the linear stability theory 
for the accelerated foil [8j, the growth time of the perturbation in the relativistic limit can 
be derived as following 



Trt m e n c A ,Xrt n 3/q / r \ 

1 b \ rrii rii L A 
Taking into account A^ ~ = 8A and L = 0.1 A in our case, we estimate that the time 
scale of the instability should be 2.2To. Such a short-wavelength perturbation grows very 
fast so that it reaches the other side of the foil soon. Finally, both interfaces are unstable and 
the entire target collapses quickly. Fig. [2](d) shows the proton energy spectrum. Although 
an energy peak is observed initially, it lowers gradually and disappears at t = 45T , leaving a 
quasi-exponential spectrum. In fact, most single-ion foils in the RPA regime show a similar 



simulation result 



141 . |15| . The main issue is the fast growth of the short- wavelength 



perturbation at the unstable interface. 

In order to check the robustness of the stable regime, we perform 3D simulations while 
keeping all the parameters same as in the 2D case except ot = 6A. A stable structure of the 
proton beam acceleration is also observed, which indicates that the regime described above 
can significantly stabilize the proton beam acceleration in the realistic three dimensional 
geometry. 

In conclusion, we present a new regime of stable proton acceleration driven by the laser 
radiation. In this regime, we smoothly extend the ID RPA model to multi-dimensional cases 
by using a two-specie ultra-thin SFT. PIC simulations show that the transverse instability 
degrades only the carbon ion acceleration and spreads them in space. The sharp front 
separating the species is always stable so that the proton layer is free from the effects of the 
RT instability. Benefiting from the superpower lasers such as HiPER and ELI, this stable 
regime might open a way to high quality proton beam generation in the future. 
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